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ABSTRACT 

We consider interacting theories with a compact internal symmetry 
group on a regular lattice. We show that the spectrum is necessarily 
vector-like provided the following conditions are satisfied: (a) weak form 
of locality, (b) relativistic continuum limit without massless bosons, and 
(c) pole-free effective vertex functions for conserved currents. 

The proof exploits the zero frequency inverse retarded propagator of 
an appropriate set of interpolating fields as an effective quadratic hamil- 
tonian, to which the Nielsen-Ninomiya theorem is applied. 



1. Introduction and Conclusions 



The only rigorous way, presently known to us, to define non-abelian gauge the- 
ories, relies on the lattice gulator. The observed fermion spectrum fits into a 
chiral representation of SU(3) xSU(2) xU(l), and so the construction of a consistent 
chiral gauge theory on the lattice has been a major goal in theoretical physics. 

In spite of extensive efforts, this program has been unsuccessful to date. The 
basic stumbhng block is the doubhng problem [1,2]. A naive discretization of the 
continuum hamiltonian of a Weyl fermion gives rise to eight Weyl fermions in the 
classical continuum limit of the lattice hamiltonian. If one starts with a Dirac fermion, 
the doublers can be eliminated by introducing the Wilson term. But the price is that 
the axial symmetry of the classical continuum hamiltonian is lost. 

There is an intimate relation between species doubling on the lattice and chiral 
anomalies in the continuum. Although classically conserved for a massless Dirac 
fermion, the conservation of the axial current is violated by quantum effects in the 
continuum theory. On the other hand, the lattice is a physical regulator, and so any 
current which is conserved at tree level should be conserved in the continuum limit 
as well. If one insists on keeping the axial symmetry of the lattice action, the conflict 
between the predictions of the lattice theory and the continuum theory is resolved 
by the appearance of the doublers. The obvious price is that now, the continuum 
limit of the lattice theory is different from the continuum theory from which we have 
started. 

Following the work of Karsten and Smit [3], the precise conditions for fermion 
doubling in a free fermionic theory defined on a regular lattice were stated by Nielsen 
and Ninomiya as a no-go theorem [4]. They assume the existence of a set of exactly 
conserved, locally defined charges which admit discrete eigenvalues. The Nielsen- 
Ninomiya (NN) theorem then asserts that there must be an equal number of positive 
helicity and negative helicity fermions in every complex representation of the symme- 
try group, provided the Fourier transform of the free hamiltonian has a continuous 
first derivative. (Recall that chirality equal helicity for massless fermions). The NN 
theorem applies in particular when the hamiltonian has a short range, and the charges 
are constructed canonically and generate a compact Lie group. 

The absence of chiral fermions is essentially a counting theorem about the zeros 
of the free hamiltonian in the Brillouin zone. It takes its simplest form in one space 
dimension. A right (left) mover is an eigenvalue of H{p) which satisfies 

E{p)^±{p-p,) + 0{{p-pcf), (1) 

in the neighbourhood of some Pc- The presence of an equal number of plus and minus 
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signs is then nothing but the famihar property that the graph of a function from the 
circle into the real numbers which has a continuous derivative and regular zeros must 
cross the axis upwards and downwards an equal number of times. 

In three space dimensions, a massless fermion corresponds to level crossing which 
is described by the effective two-by-two hamiltonian 

Heffip) = ± o-.(p - pj + 0((p - pj2) . (2) 

The point pc is called a degeneracy point. The ± signs correspond to the helicity 
of the fermion. The physically more transparent proof of the NN theorem (second 
paper of ref. [4]) relies on the examination of trajectories of constant in the energy 
surface, which go through at least one degeneracy point. Since the Brillouin zone 
is topologically a three-torus, these curves must be closed. When the curve goes 
through a degeneracy point, it crosses from positive to negative E (or vice versa). 
Moreover, one can choose an orientation to the curve using the local properties of the 
wave function, and show that near a degeneracy point the orientation (thought of as 
a unit tangent vector) coincides with the helicity of the fermion. The fact that every 
differentiable curve is orientable than implies that every crossing which corresponds 
to a left handed particle must be followed by a crossing through another degeneracy 
point which describes a right handed particle. 

The continuum limit of asymptotically free gauge theories is achieved at vanishing 
bare coupling. Together with the success of perturbative QCD in deep inelastic 
scattering, this leads to the generally accepted view that the fermionic spectrum can 
be correctly determined by setting the gauge couplings to zero. In the absence of other 
interactions, the NN theorem implies that the fermionic spectrum must be vector-like 
provided the hamiltonian has a short range. 

Attempts to avoid fermion doubling by using long range lattice derivatives lead to 
various inconsistencies at the level of weak coupling perturbation theory. Important 
examples include the SLAC derivative [5] which avoids the extra zeros by creating 
a discontinuity in the dispersion relation, and a method due to Rebbi [6] which is 
characterized by the presence of a pole in the dispersion relation. The former suffers 
from Lorentz violations and non-locality [7] while the latter suffers from the presence 
of ghosts [8] . 

As it stands, the NN theorem does not apply if the lattice model contains some 
strong interactions. This observation have led to several proposals [9] for constructing 
chiral gauge theories on the lattice which exploit a common strategy. (See ref. [10] 
for a review). One starts with a model containing only fermions and (possibly) scalar 
fields. In addition to standard quadratic terms, one introduces judiciously chosen 
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strong interactions among these fields which are operative at the lattice scale, and 
vanish in the classical continuum limit. Local symmetries of the desired continuum 
theory - the target theory - should appear at this stage as exact global symmetries 
of the model. 

The strong interactions introduced at the lattice scale do not survive in the con- 
tinuum limit. The continuum limit of the fermion-scalar model is a free theory of 
massless fermions. The hope is, that if the phase diagram is sufficiently non-trivial, 
arguments based on weak coupling perturbation theory will not apply, and a point 
in the phase diagram will be found in which the all the doublers have been decou- 
pled and the massless spectrum is chiral. It is crucial that, at the same time, the 
to-be-gauged global symmetries are not broken spontaneously. 

If this program were successful, a consistent chiral gauge theory could be ob- 
tained by turning on the gauge interactions in such a model. However, explicit model 
calculations have lead to negative conclusions in all cases studied so far [10,11]. 

Our purpose in the present paper is to provide a general treatment of the problem, 
which applies for example also to similar models based on the recently proposed 
domain wall fermions [12]. We will consider a hamiltonian defined on a regular lattice 
that has a compact global symmetry group which is not spontaneously broken. We 
will prove that under mild assumptions, which are directly related to the physical 
properties of a consistent continuum limit, the spectrum is necessarily vector- like. A 
summary of our main results is contained in ref. [13] 

Our proof is based on the observation that, while in general the physics of an 
interacting theory may be very compficated, all the crucial ingredients of the NN 
theorem have natural generalizations to this case. The object which plays the role of 
an effective hamiltonian is the zero frequency inverse propagator. Analyticity, or more 
generally, differentiability properties of the propagator are related to the range of the 
hamiltonian and hence to the locality of the theory. This relation can be exhibited 
most directly if we choose to work with the retarded propagator. Moreover, the desired 
properties of the continuum limit imply that massless fermions can still be identified 
with the zeros of the effective hamiltonian, and that these zeros still take the form of 
eq. (2). Finally, needless to say, periodicity across the Brillouin zone is mandatory in 
the interacting theory as well. 

We first present a heuristic argument that shows how, using the above ingredients, 
one can conclude that the spectrum must be vector-like in interacting theories too. 
For methodological reasons, we will base this argument on the more familiar self 
energy, i.e. the inverse of the time ordered propagator. We comment that, anyhow, 
at zero frequency all propagators coincide because the pole prescription is irrelevant. 
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For definiteness, consider a hamiltonian defined on a cubic lattice whose spectrum 
contains in particular a left handed fermion in the vicinity of the origin of the Brillouin 
zone. Suppose that this chiral fermion appears as a pole in the two point function of an 
interpolating two-component field ^. It makes no difference whether the interpolating 
field is elementary or composite. 

Consider now the self-energy of the interpolating field S~^{uj,p). According to 
our assumption, for small u and p, S^^{uj,p) ^ uj — cr-p up to higher order terms. 
The question is whether S~^{uj, p) can have no other zeros in the Brillouin zone. Let 
us consider the function cjo(p) defined implicitly by det 5'^"'^(cJo(p), p) = 0. (There 
are in fact two independent functions which satisfy this condition because is a 
two by two matrix). Notice that in general a;o(p) will be complex for real values of 
p, because in an interacting theory the self-energy has an imaginary part. 

Suppose that, For small p in the z direction, a;o(0, 0,^3) ~ p^. By the periodicity 
of the Brillouin zone, Rea;o(p) and Ima;o(p) should each have another zero as ps goes 
from zero to 27r/a. It now appears that we may be able to avoid an extra zero of 
a;o(p) by going into the complex plane. Namely, by arranging that Ima;o(p) 7^ at 
the point where Rea;o(p) — 0. This, however, is impossible! Rea;o(p) is the energy 
of the excitation, and an excitation with vanishing energy cannot have a finite width 
because of phase space arguments. Conversly, a finite width would imply that the 
vacuum is unstable! 

The above argument indicates that, as long as correlation functions have ordinary 
analytical properties, the only new possibility offered by an interacting theory is that 
of spontaneous symmetry breaking. In the present context, the scale set by the 
expectation value of the would-be Higgs field tends to be comparable to the inverse 
lattice spacing. We can try to fine tune the higgs VEV to a physical scale, but then 
the doublers will not decouple because their masses are proportional to the higgs 
VEV. For a detailed discussion of this issue see ref . [14] . 

If we forbid spontaneous symmetry breaking, we can conclude from the above 
argument that a doubler must be present. The purpose of the no-go theorem derived 
below is to cast this consideration in a rigorous form. 

We first need some definitions. Given an appropriate set of interpolating fields, 
we introduce the retarded anti- commutator 



where possible colour and flavour indices have been suppressed. We also introduce 
the space and space-time Fourier transforms 



Rap{x, t) = i9{t) (0| {V'a(x, t) , V'J(0, 0)} |0) , 



(3) 



Rap{p,t) = Y^e-'^-^Raf,{^,t), 



(4) 
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dte'^'R^f,{p,t), 



(5) 



and define 



hmRa(3{p,u; ^ ie) . 



(6) 



On the lattice one cannot demand tliat (anti)-commutators of local operators 
should vanish identically outside the light cone. We begin in sect. 2 with a discussion 
of locality properties, which are classified according to the rate at which i?Q^(x, t) 
tends to zero at large space-like separations. The main result of this section is that 
an exponentially bounded anti-commutator gives rise to an analytic 7^(p). The only 
singularities occur at generalized degeneracy points, which are those points in the 
Brillouin zone where the hamiltonian admits eigenstates of vanishing energy. The 
proof invokes the "edge of the wedge" theorem [15,16], and it is an adaptation to the 
lattice context of classic results from the theory of dispersion relations. 

In view of the intimate relation between causality and analyticity in the contin- 
uum, an exponentially bounded anti-commutator on the lattice should be a necessary 
condition for causality in the continuum limit. We comment, however, that in order to 
apply the NN theorem all that is needed is that the effective hamiltonian TZ''^{p) have 
a continuous first derivative. Consequently, it is sufficient to assume a much weaker 
form of locality, namely, that Rapi'x., t) is bounded by an appropriate inverse power of 
X. This general case will be discussed in sect. 4. Two examples of anti-commutators 
in free lattice theories are worked out in appendix A. 

In sect. 3 we characterize the possible singularities of TZ~^{p). In view of the 
desired properties of the continuum limit (before the gauge interactions are turned 
on) we may assume that at sufficiently large distances physics is correctly described by 
an effective lagrangian of massless fermions interaction only via non-renormalizable 
couplings. Demanding that all singularities of 7^~^(p) should be compatible with 
those allowed by the effective lagrangian, we show that TZ^^{p) has regular zeros of 
the form of eq. (2) which are in one-to-one correspondence with massless fermions. 
Moreover, TZ^^{p) has a continuous first derivative provided one form of singularity 
is excluded "by hand". The NN theorem is then applied to 7^~^(p) and the above 
advertized conclusions are obtained. 

The singularity we do not allow is the presence of a pole in 7^~^(p). We forbid 
this situation by assuming that the elements of the matrix TZ'^g{p) are bounded. 
We comment that, together with the assumption that symmetries are not broken 
spontaneously, this is equivalent via the Ward identities to the requirement that 
effective vertex functions, defined as correlation functions of conserved currents and 
the interpolating fields, are pole free. 
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The presence of a pole in TZ^^{p) may reflect a kinematical singularity, arising 
from a bad choice of interpolating fields. A kinematical singularity can arise, for 
example, if one uses two interpolating fields for two fermions in different corners of 
the Brillouin zone, whereas actually they can both be interpolated by a single field, 
or if not all fermions are interpolated. In appendix B we explain how one can identify 
the kinematical nature of the singularity, and describe a "trial and error" method 
for constructing an admissible set of interpolating fields which is free of kinematical 
singularities. 

If a pole in 7^^^(p) is not an artifact of an inadmissible set of interpolating fields, 
it cannot arise unless the hamiltonian is highly non-local. In a very general context, 
it has been shown that such poles give rise to the appearance of ghosts in one loop 
diagrams once gauge fields are introduced [8]. The reason is that, via the Ward 
identity, such poles appear in the vertex function, but they contribute to the vacuum 
polarization with the wrong sign. Interpreted in a hamiltonian language, this result 
implies that the action of a local current on the vacuum takes one outside the Hilbert 
space, which is unacceptable. We thus expect that it should be possible to extend our 
theorem and to rigorously exclude the presence of such poles in a consistent quantum 
theory. 

Sect. 4 is devoted to a generalization of the theorem to a much larger class of 
theories. Leaving the other assumptions unchanged, we show that a sufficient condi- 
tion for fermion doubling is that -Rq/3(x, t) be bounded by l/|x|''' where the exponent 
7 is strictly bigger than the dimension of space-time. We believe that this condition 
is also necessary, but we have not tried to construct explicit models that can verify 
this conjecture. 

In sect. 5 we discuss the implications of our no-go theorem. The theorem is suffi- 
cient to role out the existence of chiral fermions in all fermion- scalar models proposed 
so far in the literature. As already mentioned, this conclusion is in agreement with 
explicit model calculations [11]. 

Even more remarkably, the theorem implies that any attempt to reproduce the 
standard model on the lattice without violating gauge invariance must fail, if the 
spectrum can be correctly determined by switching off the Electro- Weak interactions, 
and the effective vertex functions of the Electro- Weak currents are pole-free in the 
symmetric phase. There is no need to switch off QCD in order to draw this conclusion! 
The reason is that the spectrum of QCD does not contain massless bosons, and so, 
in the absence of the photon, an effective lagrangian of the kind described above is 
valid at distances larger than one Fermi. We comment that it should be possible to 
accommodated a massless pion without changing the conclusions because a Goldstone 



7 



boson has only derivative couplings. 

The fact that gauge invariant lattice theories are necessarily vector-like raises 
an intriguing question concerning the relation between fermion doubling and the 
anomaly. If we are careful to work with an anomaly free theory and to break ex- 
plicitly at the lattice scale all global symmetries which are anomalous in the target 
continuum theory (e.g. baryon number), then there is no "need" for the appearance 
of doublers, because the lattice theory does not have a bigger symmetry compared 
to the target theory [14,17]. This is the case, for example, in the Eichten-Preskill 
model [9]. Nevertheless, the doublers do appear. We conclude this paper with a few 
comments on possible resolutions of this paradox. 

2. Locality and Analyticity 

We consider a hamiltonian Ti. defined on a three dimensional regular lattice. The 
time variable is continuous, whereas the spacial coordinates take values on a lattice 
generated by three linearly independent vectors V(fc), k = 1,2,3. We assume discrete 
translation invariance, i.e. there exist three unitary operators T(^k), commuting among 
themselves as well as with the hamiltonian, which generate finite translations 

7^fc)$(x,0^1)=$(x + V(fc),i), (7) 

where <l>(x, t) is a generic field. 

Discrete translation invariance implies the existence of a Brillouin zone B. This is 
the paralleloid spanned by the the three vectors Pq), which satisfy P(j) ■ V(fc) = 27i6jk. 
Momentum eigenstatcs arc labeled by p € and all physical quantities must be 
periodic across the Brillouin zone. 

We also assume the existence a set of conserved charges Qa which are defined 
as the sum over all lattice points of a local density. The Q-s generate the exact 
global symmetries of the model. We assume that they take discrete eigenvalues or, 
equivalently, that the global symmetry group is compact. We also assume that the 
Q-s annihilate the vacuum, i.e. there is no spontaneous symmetry breaking. 

In this paper, we will say that an operator C(x, t) is local if it depends only on the 
Heisenberg fields at time t and at sites x' whose distance from x is less than a fixed, 
finite number. Particularly relevant local operators are polynomials in the canonical 
fields, but we can in principle allow for a more general class of local operators. 

For a free lattice hamiltonian, relativistic invariance in the low energy limit 
amounts to the requirement that all fermionic eigenstates with vanishing energy com- 
pared to the lattice cutoff should be described by effective two-by-two hamiltonians 
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as in eq. (2). But if an interacting theory has massless particles in its spectrum, it 
may be difficuh (if not impossible) to separate one particle states from, say, three 
particle states containing the original particle plus a soft fermion-antifermion pair. 

On physical grounds, if the particle can be created in a causal process, there 
should exist a local interpolating field which has a finite probability to create the 
particle by acting on the vacuum. The particle should then generate a singularity in 
the two point function of the interpolating field. 

This consideration suggests the following strategy. We will not attempt to extract 
the low energy spectrum from the elusive one particle states. Instead, we will assume 
that at sufficiently large distances, physics is correctly described by an effective la- 
grangian of massless fermions interaction only via non-renormalizable couplings. We 
will then demand that all infra-red singularities of the correlation functions of the 
lattice interpolating fields should be compatible with those allowed by the effective 
lagrangian. 

We will assume that the lattice interpolating fields are local and that they belong 
to a given complex representation of the global symmetry. We intend to study their 
retarded propagator defined in eq. (3), and to establish its analytical properties. 

We first observe that i?(x, t) is bounded. This is a trivial consequence of transla- 
tion invariance and of the fact that ■0a(x, 0) is a well defined operator on the Hilbert 
space. Thus, there exists < 6i < oo such that 



Another important property is that for fixed x, i?(x, t) is an analytic function of 
t. Once more, this is a trivial consequence of the fact that on the lattice all operators 
are well defined and that ■0(x, i) = exp(i7ii)'0(x, 0) exp(— iTit). As a result, R{x.,t) 
cannot vanish identically outside the light cone, for then it would be zero everywhere. 

We should therefore distinguish between theories with various degree of locality, as 
characterized by the rate at which it!(x, t) tends to zero at large space-fike separations. 
We will say that it!(x, t) is weakly local of degree 7 > if there are positive constants 
c and 62 such that for all |x| > at 



Notice that this is an inverse power law bound. 

We will say that i?(x, t) is local if it can be bounded by an exponential, i.e. if 
there are positive constants c, 62 and /j, such that for all |x| > 



i?«^(x,t)|< V^(0,0)|0) + r{o,o)\o) <h. 



(8) 
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Finally, we will say that -R(x, t) is strongly local if it decreases faster than an expo- 
nential, i.e. if it satisfies a bound of the form (9b) for every ^. 

In this section we will assume that i?(x, t) is bounded by an exponential. This 
will allow us to prove that TZ{p) is analytic. The more general case will be treated 
in sect. 4. As suggested by the example of appendix A, it is plausible that short 
range lattice hamiltonians always give rise to strongly local (anti)-commutators. In 
any event, as we have already mentioned, an analytic propagator is compatible with 
causality in the continuum limit. We therefore expect that only theories with expo- 
nentially bounded (anti)-commutators can have a consistent continuum limit. It is 
interesting, though, that fermion doubling can be proved assuming only that -R(x, t) 
is weakly local of degree 7 > d (d is the space-time dimension). In this case Tl{p) will 
have a continuous first derivative except at generalized degeneracy points. As shown 
in appendix A, the anti-commutator of SLAC fermions violates this condition, as it 
must do in order to allow for the discontinuity of the spectrum at p = n / a. 

The last thing we need is the notion of a generalized degeneracy point. We intro- 
ducing the advanced anti-commutator 

^(x, t) = -id{-t) (0| {^,(x, t) , V^J(0, 0)} |0) , (10) 

as well as its Fourier transform A{p,ijj). For real values of p and u, we define 

Eo(p) = supjcj I ^(p, uj') = A{p, Lo') if \uj'\ < Uj} . (11) 

The physical meaning of this definition is that £^o(p) is the lowest possible energy 
for eigenstates with momentum p. We define a generalized degeneracy point by the 
condition -Eo(Pc) = 0. Thus, Pc is a generalized degeneracy point if it is the end point 
of a gap-less continuous spectrum. 

We now establish the basic analytical properties of the retarded propagator. 

Lemma. Assume that R{x,t) is local in the sense of eq. (9b). Then (a) R{p,uj) is 
holomorphic in the domain Imo; > 0, |Imp| < min{c~^Ima; , //}; (b) 7?-(p) is analytic 
with singularities only at generalized degeneracy points. 

Proof. By assumption, the r.h.s. of eq. (4) is bounded by the r.h.s. of eq. (9b) times 
g|x| |imp| fjgncg^ ^jig gyjn in gq (^4^ converges absolutely and -Ra/3(p, t) is holomorphic 
in the domain |Im p| < /x. Morever, Rapip, t) can be bounded by a polynomial of third 
degree in t times e"^*l^™Pl. The presence of the damping factor e~*^™'^ then impfies that 
the integral on the r.h.s. of eq. (5) converges absolutely for |Imp| < min{c~^Ima; , /i}. 
This proves (a). Notice that if R{x,t) is strongly local than ^^^(p, t) is an entire 
function of p. In this case R{p,u;) is holomorphic in the forward cone |Imp| < 
c~^Ima;. 
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In order to prove (b) we notice that the Fourier transform of the advanced anti- 
commutator has similar properties except that the sign of Imu; is now negative. A 
straightforward apphcation of the edge of the wedge theorem [15] now imphcs that the 
common boundary function TZ{p) is analytic, with singularities only at generalized 
degeneracy points. (A modern proof of the theorem can be found for example in 
ref. [16]). This proves (b). 

The analyticity of TZad{p) away from generalized degeneracy points, implies that 
there can be no obstructions to the smooth motion throughout the Brillouin zone 
from one zero of 7^~^(p) to another, provided we exclude the possibility of poles 
in TZ~^(p). We forbid this situation by assuming that the elements of the matrix 
7^^^(p) are bounded. The justification for this assumption has been discussed in 
the introduction. As we have explained there, we believe that this is a pathological 
situation, and that in a more complete treatment it should be possible to exclude it 
rigorously. Finally, we remind the reader that this assumption implies in particular 
the absence of kinematical singularities and hence that the set of interpolating fields 
is admissible. For more details see appendix B. 

3. A no-go theorem 

By definition, 7^~^(p) is a hermitian matrix. In order to show that it qualifies as 
an effective hamiltonian which satisfies the NN theorem, what is left for us to do is 
to show that it has a continuous first derivative at generalized degeneracy points and 
that its zeros can be identified with massless fermions. 

We assume that the continuum limit is relativistic and that the only massless par- 
ticles are fermions. Low energy physics should therefore be described by an effective 
lagrangian containing massless fermions coupled only via non-renormalizable inter- 
actions. The interaction of the smallest possible dimension is a four Fermi coupling. 
Notice that the masses of other excitations need not diverge as the lattice spacing 
tends to zero. There may well be massive excitations in the physical spectrum. In 
this case the above effective lagrangian correctly describes physics at distance scales 
which are sufficiently large compared to the Compton wave length of the least massive 
particle. 

Let us denote by Pphys the momentum variable which transforms homogeneously 
under the Lorentz group in the low energy limit. The only form of singularity in 
the zero frequency fermionic propagator which is compatible with the above effective 
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(T-Pphys (1 + Oip^pf^yjogplf^y^) 

We have Pp^jy, in front of the logarithmic term because this term involves at least 
two powers of coupling constants, and all coupling constants have a negative mass 
dimension which is at least two. 

Clearly, an allowed singularity of 7^(p) can be obtained by substituting p — Pc 
instead of Pphys is eq. (12). We call such a singularity a primary singularity. The pre- 
cise definition is as follows. A generalized degeneracy point Pc is a primary singularity 
if there exists a unitary transformation U such that 

lim [<7.(p - pj ® h] Un{p)U-' = [I^A] . (13) 

p— ►Pc 

In eq. (13), / is the identity matrix in spin space, Ii is the identity matrix in colour and 
flavour space, and ^4 is a diagonal matrix A — diag{Zi, . . . , Z^, 0, . . . , 0). The Z-s are 
non-zero constants, which are in one-to-one correspondence with massless fermions. 
The helicity of the massless fermion is determined by the sign of the corresponding 
Z. We comment that if one does not insists that the massless fermion will propagate 
at the speed of light, than one can allow for any non-singular linear transformation 
instead of U. But the relativistic case is obviously more restrictive. 

We now have to determine whether there are other allowed forms of singularity 
points. In order to do so we have to digress for a moment and discuss what are 
the acceptable values of Pc for a primary singularity point. We claim that Poincare 
invariance in the continuum limit implies that a primary singularity can occur only 
at a point which is a linear combination with rational coefficients of the three vectors 
P(j) which span the Brillouin zone. 

Assuming that this property holds and that the number of primary singularity 
points is finite (i.e. a finite number of massless fermions), there exist three smallest 
positive integers such that all primary singularity points are representatives of 
the origin of a reduced Brillouin zone B' spanned by P(j)/n(j). We can then identify 
the physical momentum with p E B'. 

Let us now consider what other forms of singularities are allowed in Tl{p). As 
an example, take a cubic lattice and assume that a primary singularity occurs at 
Pc — (27r/iVa, 0, 0) for some N. The presence of a single particle spectrum near Pc 
then implies the existence of multi-particle spectra at the points npc ior n — 2, . . . , N. 
The multi-particle thresholds will all be at £■ = because our particle is massless. 
This implies that all these points are generalized degeneracy points of the hamiltonian. 

In general, the quantum numbers of the multi-particle states may be different 
from those of the original particle. But since the global symmetry group is compact, 
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some of these states will have the same quantum numbers as the original particle and 
so they will contribute to 7?.(p) if the are no additional accidental symmetries. We will 
denote the corresponding points as secondary singularities. The leading contribution 
of the gap- less spectrum to 7?.(p) at a secondary singularity point can take the form 

± (T-Pphys Plhys Plhys ■ (14) 

As before, Pphys = P~Pc- This form is dictated by the requirement that, once we sum 
7?.(p) over all points that correspond to the same Pphys, we will obtain an expression 
compatible with the expansion of the denominator of eq. (12) in powers of the coupling 
constant. Notice that there is no reason that 7^(p) should vanish at a secondary 
singularity point, because it always receives additional, regular contributions from 
finite energy branches of the spectrum. 

In x-space, the requirement that all singularity points occur at rational points 
means that zero energy states are periodic. Eigenstates of definite physical momentum 
are then obtained by a Fourier sum over sub-lattices which respects that periodicity. 
Zero energy states are constant on every sub-lattice, and this is the lattice equiva- 
lent of the property that a relativistic continuum eigenstate of zero energy must be 
constant in space. 

By contrast, suppose that a primary singularity occurred at a non-rational point. 
This would imply the the lattice model has a non-periodic zero energy state. More- 
over, this primary singularity would lead to an infinite number of secondary singular- 
ities that would be dense on a sub-space of the Brillouin zone whose dimension is at 
least one. Under these circumstances the continuum limit cannot be both translation 
invariant and Lorentz invariant. It would be impossible to define a conserved mo- 
mentum which generates translations in space, has a continuous spectrum and which, 
at the same time, assigns the eigenvalue zero to all zero energy states. 

In fact, under these circumstances, insisting that a "momentum" operator assign 
the eigenvalue zero to all zero energy states would entail the appearance of infinitely 
many new conservation law^. Two particles whose total "momentum" is zero would 
in general never scatter into a final state with the same "momentum" and the same 
quantum numbers under the lattice internal symmetries, because in general the initial 
and final states would have different lattice momenta. Moreover, the number of 
superselection rules not explained by the lattice internal symmetries and conservation 
of this "momentum" would be infinite. 

We now collect all our intermediate results together in the following theorem. 

^ This observation is due to A. Casher. 
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Theorem. Consider a hamiltonian defined on a regular lattice. Assume the existence 
of a compact global symmetry group which is not spontaneously broken. Assume 
also that the continuum limit is relativistic and that the only massless particles are 
fermions. Under these assumptions, there is an equal number of left handed and right 
handed fermions in every complex representation of the global symmetry group, pro- 
vided i?(x, t) is local and TZ~^{p) is bounded for every admissible set of interpolating 
fields. 

Proof. Consider all sets of interpolating fields which satisfy the above assumptions 
and which belong to a given complex representation. Choose a maximal set. By this 
we mean that the total number of Z-s, as determined by the limiting procedure (13) 
and summed over all primary singularity points, is maximal. This number is then the 
total number of massless fermions in that representation. 

Locality of i?(x, t) and boundedness of TZ~^{p) imply that 7^~^(p) is analytic 
except at generalized degeneracy points. Furthermore, the allowed forms of singular- 
ities, eqs. (12) and (14), imply that 7^~^(p) has a continuous first derivative at the 
generalized degeneracy points, and that all zeros of Tl~^{p) (which occur at primary 
singularity points) are of the relativistic form (2). In addition, 7^~^(p) is hermitian. 
Hence TZ~^{p) satisfies all the assumptions of the NN theorem. Applying the theo- 
rem, we conclude that TZ~^{p) has an equal number of left handed and right handed 
zeros. Since the set of interpolating field we have chosen is maximal, this implies that 
the spectrum contains an equal number of left handed and right handed fermions in 
the given complex representation. 

For completeness, we recall why the assumption that a fermion belongs to a 
complex representation is needed in the NN theorem. If a fermion belongs to a real 
representation, it is possible to use real field (Majorana) formulation. A single Majo- 
rana fermion can then generate both a left handed and a right handed pole in its two 
point function. This is the only way to violate the one-to-one correspondence between 
poles of the two point function and massless fermions. Of course, this exceptional 
situation is of no help if we are trying to construct chiral fermions. 

4. The most general case 

We now turn to discuss the most general case. In this section we will prove that 
the spectrum is vector-like provided i?(x, t) is weakly local of degree 7 > 4, and all 
other assumptions are the same as in sect. 3. (For simplicity we continue to work in 
four dimensions). We believe that this condition is not only sufficient, but actually 
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necessary. However, we have not attempted to construct explicit models that will 
assess this conjecture. 

In complete analogy to the first part of the lemma of sect. 2, one can show that if 
i?(x, t) is weakly local of degree 7 > 3 than i?Q,^(p, t) is continuous. Moreover, since 
taking the n-th p-derivative amounts to multiplication by x" in configuration space, 
Rai3{p,t) will have continuous n-th derivatives with respect to p provided i?(x, t) is 
weakly local of degree 7 > 3 + n. Since we need a continuous first derivative, we 
assume that -R(x, t) is weakly local of degree 7 > 4. 

What is left for us to do is to prove that, if /2„^(p,t) has continuous n-th deriva- 
tives with respect to p, then the same is true for 7^(p) except at generalized degen- 
eracy points. The assumed form of the low energy effective lagrangian is the same as 
before, and so from this point on we can simply repeat the discussion of sect. 3. 

In order to relate the differential properties of /?Q^(p,t) and 7^(p) we introduce 
the spectral representation 

Rap[P,i^) = dE \ — — 

JO \ L — u E + uj 

In eq. (15) E and p are real, whereas u is complex with Imu > 0. As before, we are 
interested in the limit Ima; 0. The spectral functions are given by 

pi'j(p, E)^ [ S{E - E^)S'{p - p„) (0| ^'^(O, 0) \n) {n\ ^^{0, 0) |0) , (16) 

and p^^^(p,i?) is obtained by interchanging ipa and ipj^ and replacing 6^{p — p„) by 
6^{p + p„). In eq. (16) stands for a sum over all intermediate states. One can 
relate the two spectral functions by invoking VCT invariance, but we will not need 
this relation below. 

We will first prove the continuity of 7^(p) away from generalized degeneracy 
points. Let us denote 

^%iP)- rdEp%{p,E), i = l,2. (17) 

J u 

We want to prove the existence of a uniform bound 

I^^S(P)I<^3, (18) 

for some < 63 < cxd. Once this bound is established we are done. While we know 
very little on -Eo(p) for general values of p, our assumption that low energy physics 
is relativistic implies in particular that -Eo(p) is continuous provided its value is less 
then some physical scale Ei > 0. Therefore, for every pel?, where T> is a small 
enough neighbourhood of po € B, one has 

EmiPo) < Eo{p) , (19) 
15 




E^iPo) = min{i£;o(po) , ^^i} . (20) 

Like £'o(po), one has Em{po) = if and only if po is a generalized degeneracy point. 
Assuming that po is not a generalized degeneracy point, it is now straightforward to 
establish the following uniform bound for every p e D 

|^a/3(P,«e) - ^a/3(P,^e')l < 

^ £,''^(l^S(p.^)l+kS(p.^)l)^ 

Eq. (21) implies that TZapip) is continuous except at generalized degeneracy points. 

It remains to establish the bound (18). Since the condition of weak locahty holds 
for every t, it follows that the magnitude of long range couplings in the hamiltonian 
should obey the same bound (9a) as Ra/si'x., t), where we now set t — 0. In quantum 
mechanics, time derivatives are given by multi-commutators with the hamiltonian, 
and so the n-th time derivative of i?o,^(x, t) obeys the same bound at large space- 
like separations as Rap{'x.,t) itself. For the case at hand, 7 > 3, the x-summation 
converges uniformly for every n, and so Rap{p,t) has continuous time derivatives of 
every order. 

The spectral representation for the n-th time derivative is 



'dEE-{p'i}{p,E) + i-rp'ilip,E)} . (22) 



t=o 



Convergence of the integral on the r.h.s. of eq. (22) for even n imphes that \Pap{p, E) + 
p^^p(p,E)\ decreases faster than any power of £■ for £■ — > 00. Convergence of the 
integral for odd n implies the same for the difference of the spectral functions. Hence 
each of the spectral functions separately decreases faster than any power. Together 
with continuity of Rap{p, t) as a, function of p, this implies that Q^*^(p) as defined in 
eq. (17) is continuous. Compactness of the Brillouin zone than implies the existence 
of the uniform bound (18). This concludes the proof that TZ{p) is continuous except 
at generalized degeneracy points. 

The proof that the derivative of TZ{p) is continuous except at generalized de- 
generacy points is completely analogous. We simply replace every function of p by 
its p-derivative in eqs. (15-22). In general, this process can continue up to the n-th 
p-derivative, provided i?a/3(x, t) is weakly local of degree 7 > 3 + n. 

In conclusion, we showed that if Rapi'^jt) is weakly local of degree 7 > 4 than 
TZ{p) has a continuous first derivative except at generalized degeneracy points. From 
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here on one can repeat the discussion of sect. 3 and prove that the spectrum is 
necessarily vector-hke. 

5. Discussion 

The main lesson of our no-go theorem is that, as far as fermion doubling is con- 
cerned, there is essentially no difference between free and interacting lattice theories. 
From the correlation functions of the interacting theory one can always construct 
an object which plays the role of an effective quadratic hamiltonian for the massless 
fermions. The types of singularities in the dispersion relation which allow one to 
escape the conclusions of the no-go theorem are the same as in the free case. Thus, 
the only remaining question is whether such "bad" singularities can be less danger- 
ous if they occur in an interacting theory. In this paper wc have gone a long way 
towards giving a rigorous, negative answer to this question. For example, we can 
reject without relying on perturbation theory [18] any lattice model whose dispersion 
relation for the effective low energy degrees of freedom contains a discontinuity, on 
the grounds that such a theory must be a-causal. 

Our no-go theorem is sufficiently strong to exclude the existence of chiral fermions 
in all the fermion- scalar models proposed for that purpose in the literature [9]. This 
is no surprise, as all explicit model calculations [11] have consistently reached the 
same conclusion. The virtue of our approach is that it provides a uniform basis 
for treating all these models, and that it constraints the spectrum by invoking only 
physical properties which characterizes a consistent continuum limit. 

As we have already mentioned in the introduction, perhaps the most striking 
consequence of our theorem is the constraints it puts on any attempt to reproduce 
the standard model on the lattice without violating gauge invariance. It asserts that 
any such attempt must fail, if the spectrum can be correctly determined by switching 
off the Electro- Weak interactions (there is no need to switch off QCD), and provided 
that the effective vertex functions of the Electro- Weak currents are pole-free in the 
symmetric phase. 

One may question the relevance of this conclusion on the grounds that we are not 
dealing directly with gauged Electro- Weak interactions. But when one says that the 
a gauge theory is chiral, one is making a statement on its elementary, gauge variant 
fields! These cannot be related directly to physical observables without invoking 
additional assumptions such as complementarity. Thus, the very definition of a chiral 
gauge theory implicitly assumes that we can determine the elementary fermions' 
content by setting the gauge coupling to zero. 
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Moreover, the widely accepted view that QCD, as defined on the lattice, is a 
consistent quantum theory with a non-trivial continuum limit, is based on asymptotic 
freedom and the validity of weak coupling perturbation theory at short distances. 
Thus, in our opinion, one cannot reject our conclusions without, at the same time, 
questioning the validity of the present understanding of non-abelian gauge theories. 

From a conceptual point of view, this situation is very intriguing, for Electro- Weak 
phenomenology has been found to be in remarkable agreement with the predictions 
of continuum perturbation theory, and the lattice has been introduced mainly to 
allow for a quantitative treatment of the strong interactions. Here we are faced with 
an inability to provide a uniform, consistent treatment to all the interactions of the 
standard model simultaneously. 

Paradoxically, the only approach to putting the weak interactions on the lattice 
which still seems viable is based on using gauge non-invariant actions. Several such 
proposals exist in the literature [19]. We include among them the Rome approach [20], 
because one cannot derive the gauge fixed action of this approach from a gauge 
invariant lattice action. It has also been claimed [21] that the recently proposed 
domain wall fermions do not give rise to a chiral spectrum unless one gives up tree 
level gauge invariance. 

So far, there are only partial results concerning these models, and all of them are 
in the context of perturbation theory. At this level, it is not unreasonable to expect 
that gauge invariance will be recovered in the continuum limit, provided the spectrum 
is anomaly free. In this respect the lattice is not very different from continuum 
field theories, where no gauge invariant regularization for chiral theories is known 
to exist. Nevertheless, the Adler-Bardeen theorem [22] gauranties the consistency of 
perturbation theory provided the theory is anomaly free at the one loop level. In 
modern language, the Adler-Bardeen theorem can be derived using renormalization 
group arguments [23]. The renormalization group exists on the lattice as well, and 
so it is not unreasonable to expect that some version of the Adler-Bardeen theorem 
should hold on the lattice. 

The real unknown is the behaviour of such models at the non-perturbative level. 
By the no-go theorems, gauge invariance on the lattice implies a vector-like spectrum 
if we demand that the continuum limit be free of inconsistencies such as violations 
of causality or Lorentz invariance. Thus, one can suspect that models with explicitly 
broken gauge symmetries will reveal some unexpected features at the non-perturbative 
level. Until a detailed understanding of non-perturbative effects in such models is 
reached, their consistency as well as their relevance to the real world remain unclear. 
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Appendix A 

In this appendix we give two examples of anti-commutator functions in free 
fermionic hamiltonians in one space dimension. The retarded anti-commutator is 
obtained by simply multiplying the ordinary anti-commutator by 0{t). 

Let us define 

Af{x, t) = {0\mx, t) , z^t(o, 0)}|0) . {A.l) 

Here ijj{x.t) is a single component fermion field. In this appendix and the next one 
we set a = 1, and the space coordinate x takes integer values. AF{x,t) satisfies the 
homogeneous wave equation 

{dt-V)AF{x,t)^0, (A2) 

where V is a lattice difference operator, and the boundary condition Af{x, 0) = S^fi- 
It is convenient to express Aj?(x, t) as 

AF{x,t) = -i{dt + V)A{x,t), (A3) 

where A{x, t) is a homogeneous solution of the second order equation 

(d^-V')A{x,t)^0, (A4) 

with the boundary conditions A{x, 0) = and dtA{x, 0) = i6x,o- An explicit repre- 
sentation for A{x,t) is 

A{x,t) = ^ [dcu r dpe^^-'P'^Siuj^ - E^ip))e{uj) (A.5a) 

277 J-n2\E{p)\ V J ^ ' 

Here E{p) is the dispersion relation as determined by the free hamiltonian. 

Consider first a nearest neighbour hamiltonian. The lattice difference operator is 
V = ^{Sx,x'+i — (^x,x'-i), and the dispersion relation is E{p) = sinp. It is convenient 
to calculate dtA{x,t) first. A straightforward calculation gives rise to 

d A( t) { ^ ^ X = 2n -|- 1 {A 6) 

iJ\x\{t) , x^2n. 
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Using eq. (A. 3) and the boundary conditions satisfied by A{x, t) we obtain 

A,(x,t) = | fl^l(')' ^ (A7) 

This result clearly exhibits the presence of a primary right mover and a doubler which 
is a left mover. We can estimate AF{x,t) using the asymptotic expansion for Bessel 
functions where both the order \x\ and the argument t are large. Denote ^ = {\x\—t)/t. 
Assuming ^ <C 1 and ^ 1 one has 



J\x\{t) 



cos 



The slow decay of Ap{x,t) inside the light cone signals the presence the the high 
energy excitations with group velocity which is smaller than one. Outside the light 
cone, Aj?(a;, t) decreases faster than an exponential. 

We now perform a similar calculation for the SLAG derivative. The dispersion 
relation is E{p) = p. We now have 



Af{x, t) = ^ fdu r dp {cu + p)e*'"*-^^'^5(u;2 - E'^ip))e{uj) 

ATT J J-TT 

1 sin(7r(x — t)) 



TT X — t 



(A.9) 



Thus, although A{x,t) is maximal on the classical trajectory x — t, it decays ex- 
tremely slowly for both \x\ < t and |x| > t. The non-local character of the anti- 
commutator is necessary in order to produce the expected discontinuity in TZ(p). 



Appendix B 

An unwise choice of interpolating fields may give rise to two types of spurious, 
kinematical singularities in 7^^^(p). This can best be illustrated through an example 
in the context of a free fermion theory. Consider again a single component field in 
one space dimension with a standard nearest neighbour hamiltonian. The spectrum 
consists of a right mover at p = 0, and a left mover at p = tt. In self-explanatory 
notation, the retarded anti-commutator is 

7^(p) = {^^^){p) 
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Suppose that, instead of using i(j{x) as a single interpolating field for both fermions, 
we unwisely decide to use the two fields ^±{x) = ±^{'ijj{x + 1) + "ipi^x — 1) ± 2iIj{x)). 
In Fourier space this becomes ip±{p) = (1 ± cosp)iJj{p). The ± is now considered as 
a "flavour" index, and 



R'(x, t) = 
In Fourier space this becomes 



{B.2) 



[1 + cosp)^ sin^p 

sin^p (1-cosp)^ / sinp 



Notice that is nothing but TZ{p). We now observe that det7?.'(p) = identically. 
This is a consequence of the fact the the two fields ip±{x) are linearly dependent in 
the Hilbert space. By this we mean that their matrix elements between any two states 
of momentum k and k + p are linearly dependent for fixed p. In the above example, 
the proportionality constant is the kincmatical factor (1 + cosp)/(l — cosp). 

The solution to this problem is simple. In the above example, we have to replace 
the two fields tp±{x) by a linear combination of them aip+{x) + htp-{x). In so doing, 
there are some constraints on the coefficients a and b which are necessary in order to 
prevent the appearance of spurious poles in 7^(p). 

Clearly, such a pole will appear if one of the coefficients is zero. In this case we 
are interpolating only one of the fermions, and a zero in the propagator will appear 
at the location of the missing fermion. Another possibility is that we take both a and 
b to be non-zero, but with opposite signs. What is common to both is the existence 
of a point where the linear combination a(l + cosp) + 6(1 — cosp) vanishes. On the 
other hand, any choice of a and b where both are, say, strictly positive, is admissible. 
In this case the kinematical factor never vanishes and both fermions are interpolated. 
Notice that the choice a = 6 = 1/2 brings us back to the original field ■ip{x). 

Of course, similar phenomena can arise from a bad choice of interpolating fields 
in an interacting theory as well. What we need is a procedure that will allow us 
eliminate all kinematical singularities. Since we have no knowledge on the way a 
given set of interpolating fields was constructed, we are lead to the following two step 
procedure. 

In the first step we on purpose multiply the interpolating fields by various kine- 
matical factors. If necessary, we enlarged the set by additional fields that are linearly 
dependent in the Hilbert space on existing ones. The aim at this stage is to achieve 
a one-to-one correspondence between fields and fermions. 
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In the second stage we minimize the total number of interpolating fields by taking 
appropriate linear combinations of the previous, linearly dependent fields. Suppose 
that there are n primary singularities, and A:„ fermions are interpolated at the n-ih. 
primary singularity. We can than replace the fields which interpolate the first fermion 
at every point by a single linear combination. This single field will interpolate the 
first fermion at every primary singularity. 

If we continue this process, we are guaranteed that the final set will be inde- 
pendent in the Hilbert space. However, without additional information, we do not 
know whether in every step it will be possible to avoid the appearance of zeros in 
the propagator. If no choice of interpolating field can be made which is free of such 
zeros, than these zeros are genuine, and the pole in the inverse propagator is not a 
kinematical singularity. As we have discussed in the introduction, we expect that 
theories with this property must be inconsistent. 

Finally, we comment that the complicated process described above was necessary 
only because we assumed that we know nothing about the hamiltonian of the theory, 
nor on the way the interpolating fields are constructed from the elementary fields 
of the theory. In practise, one has a specific model in mind, and it is easy to iden- 
tify reasonable candidate interpolating fields as well as to verify that no kinematical 
singularities are present. 
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